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^ . Abstract 
Or 

^ ^ We compare different analytical and numerical methods for studying the par- 

> ■ titions of a finite system into fragments. We propose a new numerical method of 

exploring the partition space by generating the Markov chains of partitions based 
on the Metropolis algorithm. The advantages of the new method for the problems 
where partitions are sampled with non-trivial weights are demonstrated. 



PACS numbers: 25.70.Pq, 02.70.Lq . 

Many fields of physics deal with the common phenomenon that, under appropriate 
conditions, a compound system can disintegrate into constituents. Let us consider an 
isolated system composed of Aq identical particles (we call them nucleons) which are kept 
together by some attractive forces. If sufficient energy is put into the system, it will 
disintegrate into fragments. These fragments can either be individual nucleons or bound 
clusters of several nucleons. Examples of such processes abound in condensed matter 
physics, nuclear physics, and astrophysics. In order to provide a microscopic description 
of such processes, one must sort out possible partitions of the system and compare their 
probabilities. At the first step, it is necessary to develop methods of generating and 



sampling the partitions. The aim of this paper is to propose a new and efficient method 
of doing this. 

The obvious way to proceed is simply to construct all partitions directly and calculate 
the characteristics of interest. Unfortunately, this approach is possible only for small 
Aq because the total number of partitions, P(Aq), grows rapidly with Aq. For instance, 
P(100) = 190569292 while P(200) = 3972999029388. Even if one needs only perform a 
few non-trivial operations for each partition, this task becomes intractable for Aq > 100. 
We shall, however, reserve this direct method for checking the more practical methods 
presented below. 

First, we address an analytical approach to dealing with the Euler's partitioning prob- 
lem. It is based on the Generating Function (GF) formalism ||T|. This approach can be 
applied successfully for calculating average characteristics of partitions. We characterize 
each partition / by the multiplicities {Na} of fragments with different nucleon numbers A, 
1 < A < A . Then, the conservation of the total nucleon number for each / is expressed 
as: 

A 

E N A ] A = Aq . (1) 
A=l 

Evidently, the total fragment multiplicity M in the channel / is 

A 



M f = Yl N A f) . (2) 
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Following a well-established method in mathematical literature 0], we introduce an un- 
constrained generating function Z(x): 
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where the ca are arbitrary numbers which can later be taken as 0^=1. Here x can be 
considered as a Lagrange multiplier. Now we can calculate the total number of partitions, 
P(Aq), by simply expanding eq.(Q) and counting the coefficient of x A °. The results for 
large Aq or x —>■ 1 are well approximated by famous Hardy-Ramanujan formula: 
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One can use this generating function to calculate approximately the average multiplicities 
of fragments (Na) over all partitions. This is done by replacing the exact constraint of 
eq. (P by an approximate one: 

oo 

J2 (N A )A = A , (5) 

A=l 

i.e. the constraint is fulfilled on average only. Then one obtains 

dhx(Z(x)) ~ Ax A 



Aq = x- 



dx v , 



S —a . (6) 
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where we have set the Ca — 1- This equation must be solved to determine x. A very good 
approximation to the solution at large A is 



Now, the mean multiplicities of fragments can be calculated as: 

The result is shown in fig. 1 (top panel) in comparison with the results of the direct 
method in which all the partitions are included in the calculation. It is seen that the 
agreement is good except for a slight discrepancy at large A which indicates an expected 
finite size effect. Indeed, eq. (§) gives small but finite (Na) even for A > A when the 
exact calculation gives strictly zero. The average multiplicity of all fragments can be 
calculated as (M) = J2a(Na) and is well approximated by the expression 



with b =0.315087. For example, for y4 =100 it gives us (M) = 21.32 while the exact value 
obtained with the direct method is 21.75. 

More generally, it is useful to consider the situation in which partitions are biased 
with certain weights. In statistical theory, for example, identical fragments are counted 
in a partition sum with a factorial weight l/iV^!. The weight of a partition is then 
Wf = 1/ 11a Na- In this case, the corresponding generating function can be written as: 
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This form is similar to the grand canonical partition sum if one identifies x with the 
fugacity and ca's with the internal partition sums of individual fragments [||. Now instead 
of eqs.djD and (||) one easily obtains (after substituting 0^=1): 

oo 

A =J2Ax A , (N A )=x A . (11) 

A=l 

For Aq — > oo one finds the approximate expressions x = exp(— I/^/Aq) and (M) = \TA^. 
These results are shown in fig. 1 (bottom panel). The mean multiplicity (M) = 10 for 
the case A Q = 100 is in good agreement with the exact value of 9.77 obtained by direct 
calculation. 

For the two simple examples considered above one can calculate also the multiplicity 
distributions of individual fragments. It is clear from the structure of the generating 
functions, eqs. (f|) and fllQl) , that the distribution is exponential in the first case and 
Poissonian in the second case. The normalized multiplicity distributions are respectively, 

As seen in fig. 3, the exact results are reproduced by these distributions with high 
accuracy. 

In practice, however, direct accounting for all partitions can only be done for A < 100. 
If the weight factors are complicated, it can also be hard to find an analytical solution. 
Multiplicity distributions and correlations, which are of considerable physical interest, are 
particularly difficult to obtainP]. There is thus a need for another method, presumably 
based on the generation of individual partitions. Obviously, it must be efficient enough 
to permit computer simulation within a reasonable time. 

A first attempt to develop such a method was made in refs. 0] by introducing a bias 
function b(A , M) = P(A , M)/P(A ), where P(A , M) is the total number of partitions 
with exactly M fragments. It can be calculated using the recursion relation |], | 

P(A , M) = P(A — M, M) + P(A — 1, M — 1) . (13) 

this respect an interesting development of an analytical method was recently made in ref. [[?). 
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As before, the total number of partitions is 

P(4))=£P(A ,M). (14) 

M 

This bias function is used to generate a sample of partitions by the Monte Carlo method. 
First, M is selected randomly with a probability given by the bias function, b(Ao,M). 
Then, a random partition with selected multiplicity is generated as described in ref. 0. 
We shall refer to this method as Biased Random Generation (BRG) . Another Monte Carlo 
method of generating partition samples using a bias function obtained with a Laplace 
transformation is described in ref. ||. 

Figs. 1 and 2 (top panel) show how well the BRG method works in the case when all 
partitions have equal weighs. The results are presented for t4o=100 and summarize the 
outcome of 10 5 randomly generated partitions. By construction, this method is guaranteed 
to give the correct multiplicity distribution as shown in fig. 2 (top panel). It is less 
trivial that it reproduces correctly also the mean multiplicities of individual fragments 
as well as other distributions. Unfortunately, the BRG method has a serious drawback: 
It produces correct results only for the case in which the weights of partitions are equal. 
This is not surprising given that eq. (|l^) was obtained under this assumption. When we 
introduce nontrivial weight factors, for instance relative factorial weights W = 1/ILi-Wa! 
for partitions with fixed M, the method fails. This is clearly seen in the bottom panel 
of fig. 1 for mean fragment multiplicities. In the case of nontrivial partition weights 
the calculation of a bias function might be even more difficult than the calculation of a 
corresponding generating function. 

Here, we propose a new method of the partition sampling which is designed especially 
for computer simulations. The idea is to generate a Markov chain by moving from one 
partition to another by minimal steps, i.e., by demanding that neighboring partitions differ 
by the state of one nucleon only. We shall refer to this method of generating partition 
samples as Markov Chain Generation (MCG). The procedure allows the following moves: 
(a) to transfer a nucleon from one fragment to another, (b) to make a nucleon free, or 
(c) to attach a free nucleon to a fragment. In addition, one must ensure that each new 
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partition is different from the previous one, since fragments with the same A are to be 
regarded as indistinguishable. 

As well known, any sampling procedure of this kind must satisfy the detailed balance 
requirement. This can be achieved by applying the famous Metropolis algorithm |J, 
where a chain of partitions is generated by performing subsequent moves in the partition 
space biased by the partition probabilities W (weight factors). As shown elsewhere (e.g., 
ref. 0), this method provides a correct description of the complete partition space for 
any specified weight factors W. In the MCG the number of all possible moves is limited 
and easily countable for any partition. By generating a new partition we account for the 
probability of all possible moves, and thus we avoid the bias function problem. Detailed 
balance is guaranteed by application of the Metropolis algorithm. 

The numerical procedure is implemented in the following way: 
Step I: For a given partition with M fragments of mass numbers A; (i = 1, ...,M), enu- 
merate all fragments in the order of decreasing mass so that A\ > A 2 > ... > Am- This 
order is to be strictly maintained; any move violating this ordering is rejected. In this 
manner, we ensure that each move gives a genuinely new partition. 
Step II: Select at random the fragment % that looses a nucleon and the fragment j 
(j = 1,...,M + 1; j ^ i) that accepts it. (The case j — M + 1 corresponds to mak- 
ing the nucleon free.) Check this move against the ordering requirement of Step I. If the 
order is violated, repeat the determination of i and j. 

Step III: Calculate the weight of a new partition, W new , and compare it with the weight 
of the previous one, W Q id- A new partition is added to the ensemble if W new > W \d- If 
W new < W Q \d, a new partition is added with probability W new /W id- Otherwise, the old 
partition is taken as the new one and a new move is undertaken. 

Step IV: Calculate the characteristics of interest by taking all partitions from the chain. 
The chain is truncated when these characteristics are saturated. 

We stress that, contrary to the GF and BRG methods discussed above, the MCG 
method is a purely numerical procedure which requires nothing more than random number 
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generation. This provides a welcome degree of universality which is missing in other 
methods. For example, similar to the direct calculation, our method can be applied in 
case of any partition weight, as well as it can be easily generalized for other partition 
spaces, e.g., when fragments are characterized by two numbers (such as mass A and 
charge Z) instead of one ||. 

The initialization problem, i.e., the question of which partition should be taken as a 
seed, does not appear to be important for the MCG method. The system with Ao=100 
loses all memory of the initial partition after approximately 10 4 moves. In order to obtain 
a representative partition sample, one should just discard these initial partitions from the 
ensemble. This is verified for several cases when partition weights vary smoothly with 
fragment mass and the number of fragments. In other cases, the number of initial moves 
may increase. This problem must be analyzed in each particular case. 

We have checked the MCG method in a number of ways. The results are presented in 
figs. 1-4 for two cases: first, when all partitions have equal weights (top panels) and second, 
when partitions with identical fragments are suppressed by the factorial weights Wf = 
1 / JIa (bottom panels). They show the mean fragment multiplicity as a function of A 
(the mass distribution), the distribution of total fragment multiplicity, and a very specific 
characteristic, i.e., the distribution of multiplicities of particular fragments (A = 1, A = 4 
and A > 10) taken over all partitions. The results of the exact direct method and of 
Markov chain generation are in remarkably good agreement. It should be stressed that 
for Aq = 100 all 1.9 x 10 8 partitions are included in the direct method while only 10 5 
partitions can be taken from the chain to explore the entire partition space with the 
MCG method. Small discrepancies in the tails of the distributions are seemingly related 
to a limited sample size and numerical precision. However they are not important in 
practice because of their very small relative weight in the chain. For smaller systems 
(e.g. A =20) the agreement is also good. For larger systems (e.g., A =1000), where the 
direct method is intractable, comparisons were made with the analytical GF method. As 
demonstrated in fig. 4, the agreement is quite good, apart of a small discrepancy in the 
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tails. One should bear in mind, however, that the GF method slightly overestimates the 
exact result (see fig. 1). We emphasize that the same high quality agreement between the 
direct and MCG methods is achieved in both considered cases which differ significantly by 
the weight factors. Calculations have been made for partitioning with other weights, and 
similar agreement has been found. Therefore, we believe that the MCG method described 
here offers a simple and correct numerical solution to the partition sampling problem. 

In conclusion, we have analyzed several methods for calculating characteristics of the 
partition space of a finite composite system. We have developed a new numerical method, 
the Markov Chain Generation, which is flexible and efficient in practical calculations 
with complicated partition weights. We see a variety of applications of this method in 
different fields dealing with finite-size objects, from atomic nuclei to molecular clusters 
and astrophysical objects. We believe that this method will be very useful for studying 
the thermodynamics of finite systems p], ||, |5|, f§. 
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Figure captions 



Fig. 1. Average multiplicities (N A ) of fragments with mass number A for the system 
with total mass A = 100. Solid lines: direct calculation taking into account all parti- 
tions, dashed lines: numerical Markov chain generation of partitions, dot-dashed lines: 
analytical calculations by the generating function method, dotted lines: biased random 
generation. Top panel: for partitions with equal weights, bottom panel: for partitions 
with the factorial weights 1/ U A N A \. 

Fig. 2. Distribution of total fragment multiplicity M for the system A = 100. 
Notations are the same as in fig. 1. 

Fig. 3. Multiplicity distributions of fragments with A = 1, A = 4 and A > 10 for the 
system A = 100. Notations are the same as in fig. 1. 

Fig. 4. Comparison of fragment mass distributions for A =20 and 1000 calculated 
by the analytical and Markov chain generation methods. Top and bottom panels show 
calculations for two different weighting factors as above. 
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